Metallic surface electronic state in half-Heusler compounds iiPtBi (R = Lu, Dy, Gd) 
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Rare-earth platinum bismuth (i?PtBi) has been recently proposed to be a potential topological 
insulator. In this paper we present measurements of the metallic surface electronic structure in 
three members of this family, using angle resolved photoemission spectroscopy (ARPES). Our data 
shows clear spin-orbit splitting of the surface bands and the Kramers' degeneracy of spins at the F 
and M points, which is nicely reproduced with our full-potential augmented plane wave calculation 
for a surface electronic state. No direct indication of topologically non-trivial behavior is detected, 
except for a weak Fermi crossing detected in close vicinity to the F point, making the total number of 
Fermi crossings odd. In the surface band calculation, however, this crossing is explained by another 
Kramers' pair where the two splitting bands are very close to each other. The classification of this 
family of materials as topological insulators remains an open question. 



I. INTRODUCTION 

The discovery of topologically non-trivial states of 
matter opens up a new realm of knowledge for funda- 
mental condensed matter physics. Unlike conventional 
materials, these "topological insulators" exhibit metallic 
surface states that are protected by time reversal sym- 
metry, while maintaining an insulating bulk electronic 
structure. This leads to a variety of novel properties in- 
cluding odd number of surface Dirac fermions, strict pro- 
hibition of back-scattering, etc., paving the way to poten- 
tial technical breakthroughs in e.g. quantum computing 
process via the application of spintronics^ ^. Recently, 
extensive theoretical and experimental efforts have led 
to the realization of such fascinating behaviors in e.g. 
the HgTe quantum wells^^"^, the Bii-^^Sb^; system^^^and 
the Bi2X3 (X = Te, Se) binary compounds^ Numer- 
ous half-Heusler ternary compounds have been proposed, 
theoretically, to be pote ntial n ew platforms for topolog- 
ical quantum phenomenally^, where the inherent flexi- 
bility of crystallographic, electronic and superconducting 
parameters provide a multidimensional basis for both sci- 
entific and technical exploration. The experimental de- 
termination of their topological class would set the basis 
for possible spintronic utilization and further studies on 
the interplay between the topological quantum phenom- 
ena versus e.g. the magnetic^^, superconducting^^ and 
heavy FermioniJI^ behaviors. 

Theoretically, the topological insulators experience a 
gapless surface state protected by time reversal symme- 
try and thus are robust against scattering from local im- 
purities. Such a surface state is "one half of a nor- 
mal metal in that the surface bands are strongly spin- 
polarized, forming a unique spin helical texture^ On 
the other hand, the Kramers' theorem requires that the 
spin be degenerate at the Kramers' points - /c-points of 
the surface Brillouin zone where time reversal symmetry 
is preservecP^. At the interface between, say, a normal 
spin-orbit system and vacuum, the spin-polarized sur- 



face bands connect pairwise (Kramers' pair), crossing the 
chemical potential /i an even number of times between 
two distinct Kramers' points. At the interface between a 
topologically non-trivial material and vacuum, however, 
one expects the surface bands to cross /i an odd number 
of times^. 

In this paper we present a systematic survey on the 
surface electronic structure of half-Heusler compounds 
i?PtBi {R = Lu, Dy, Gd) using angle resolved photoe- 
mission spectroscopy (ARPES). Our results show clear 
spin-orbit splitting of the surface bands that cross the 
chemical potential, which is nicely reproduced in the full- 
potential augmented plane wave calculation for a sur- 
face electronic state. The Kramers' degeneracy of spin 
is unambiguously detected at both the F and M points. 
No direct indication of topologically non-trivial behav- 
ior is detected, except for the fact that there is a weak 
Fermi crossing in the close vicinity to the F point, making 
the total number of crossings five. In the surface band 
calculation, however, this inner crossing is explained by 
two spin-orbit splitting bands that are very close to each 
other, forming another Kramers' pair. In this band con- 
figuration, the total Berry phase would be zero for the 
half-Heusler systems, and they would not be topologi- 
cally non-trivial. The detailed topological class of this 
family of materials thus remains an open question, re- 
quiring a detailed spin-resolved ARPES study with ultra- 
high momentum resolution and a direct calculation of the 
topological invariants based on the first principle band 
structure. 



II. EXPERIMENTAL 

Single crystals of KPiBi {R = Lu, Dy, Gd) were grown 
out of a Bi flux and characterized by r oom temperature 
power X-ray diffraction measurementS"*^^'^^. The crys- 
tals grow as partial octahedra with the (111) facets ex- 
posed. Typical dimensions of a single crystal are about 
0.5 X 0.5 X 0.5 mm^. The ARPES measurements were 
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FIG. 1: (Color online) Surface Fermi maps of half-Heusler compounds i?PtBi {R = Lu, Dy, Gd). (a) Clb Crystal structure 
of i?PtBi. The crystallographic axes are rotated so that the (111) direction points along z. The red parallelogram marks the 
Bi(lll) cleaving plane, (b) The surface and bulk Brillouin zone for the rotated crystal structure in (a). Here kz corresponds 
to the (111) direction of the fee Brillouin zone, (c)-(e) Surface Fermi maps of i?PtBi. All data is taken with 48 eV photons at 
T — 15 K. Yellow lines denote the surface Brillouin zone. 
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FIG. 2: (Color online) Surface electronic structure of GdPtBi: 
Comparison between ARPES data and calculational result, 
(a) Fermi map of GdPtBi observed by ARPES, same as Fig. 
1(e). (b) Calculational surface Fermi map of GdPtBi at the 
Bi(lll) cleaving plane. See text for details, (c) ARPES band 
structure along the contour F-M-^-F. Inset of (c) enhanced 
the ARPES intensity near M and K for better visibility of 
the bands, (d) Calculational band structure with respect to 
(c). (e)-(f) Expanded figures for (b) and (d), respectively, 
showing six Fermi crossings. Panel (e) is rotated by 30° with 
respect to (b). 



performed at beamline 10.0.1 of the Advanced Light 
Source (ALS), Berkeley, California using a Scienta R4000 
electron analyzer. Vacuum conditions were better than 



3 X 10"^^ torr. All ARPES data was taken at T = 
15 K, above the magnetic ordering temperatures of all 
compound The energy resolution was set at ~ 15 
meV. All samples were cleaved in situ, yielding clean 
(111) surfaces in which atoms arrange in a hexagonal 
lattice. High symmetry points for the surface Brillouin 
zone are defined as r(0,0), K{ko,0) and M(0,/coV3/2) 
with unit momentum ko = V^Tr/a, where a is the lat- 
tice constant for each type of crystals. We emphasize 
here that no stress or pulling force is felt by the sam- 
ples, which ensures that the measured data reveals the 
intrinsic electronic structure of the single crystals. 



III. RESULTS AND DISCUSSION 

We begin this survey in Fig. 1 by showing the Fermi 
maps of the three half-Heusler compounds i?PtBi {R 
= Lu, Dy, Gd). Previous theoret ical calculations for 
the bulk electronic structur^^^^^^I^ suggested that the 
Kramers' crossing at the f point happens very close to /i; 
the Fermi surface reduces to a single point (Dirac point) 
at f . The data in Fig. 1 shows that, at least in the 
(111) cleaving plane, this is not the case. Instead there 
are several bands crossing /i in the vicinity of both the 
f and M points. The overall Fermi surface for all three 
half-Heusler compounds are similar, indicating a similar 
cleaving plane and band structure for all members. By 
comparing the band structure measured at the (111) sur- 
face with results of band calculations for GdPtBi (Fig. 
2), we find the cleaving plane to be Bi(lll), marked by 
a red parallelogram in Fig. 1(a). A closer look at Fig. 
1(c)- (e) reveals that the f pockets have different sizes 
for different Heusler members. For example the circu- 
lar F pockets in LuPtBi are larger in size than those in 
GdPtBi. This indicates a different effective electron oc- 
cupancy for different members of the half-Heusler fam- 
ily. One should also note that in Fig. 1(e) the inner of 
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FIG. 3: (Color online) Absence of kz disper- 
sion as proof for the observation of a surface 
electronic structure, (a)-(b) kz (r-A) disper- 
sion maps for LuPtBi. Data is obtained by 
scanning the incident photon energy hu from 
30 to 80 eV along (a), the T-K and (b) the f- 
M direction, (c) Calculational Fermi surface 
map for the bulk state in the A-T-K plane. 
See panel (a) for comparison, (d)-(g) Band 
dispersion maps along the the T-K direction 
for selected hus. It is clear that all observed 
bands are independent of hu (kz). (h) De- 
tailed peak analysis for the momentum dis- 
tribution curves (MDCs) at the chemical po- 
tential for four different photon energies, k 
range is indicated by an orange double arrow 
in (d). Bars in different colors indicate the 
Fermi crossings for different bands. 



the two bright f pockets is hexagonal in shape, reminis- 
cent of the hexagonal shape of the Dirac cone in Bi2Te3 
(Ref.^), which is explained by higher order terms in the 
k ' p Hamiltonian^^. This hexagonal shape is very nicely 
reproduced in the calculation [Fig. 2(b)]. For clarifying 
the topological class of the half-Heuslers, two immediate 
questions follow the observations in Fig. 1: (1) Are the 
observed bands actually arising due to the sample sur- 
face? (2) Exactly how many times do the bands intersect 
the chemical potential along the f-M line segment? 

Fig. 2 shows the comparison between the ARPES data 
and a calculational surface state in GdPtBi. For both the 
band structure and Fermi surface calculation, we used a 
full-potential linear augmented plane wave (FPLAPW) 
method^^ with a local density functional^^. The crys- 
tallographic unit cell is generated such that the (111) 
direction of the fee Brillouin zone points along the z- 
axis. For calculation of the surface electronic structure, 
supercells with three unit cell layers and a 21.87 a.u. vac- 
uum is constructed. We calculated band structures of all 
six possible surface endings (Gd-Bi-Pt-bulk, Gd-Pt-Bi- 
bulk, Bi-Gd-Pt-bulk, Bi-Pt-Gd-bulk, Pt-Gd-Bi-bulk, Pt- 
Bi-Gd-bulk); only the Bi-Pt-Gd-bulk construction shows 
good agreement with experiment [Fig. 2(b), (d)-(f)]. 
Structural data were taken from a reported experimen- 
tal result^^. To obtain the self-consistent charge density, 
we chose 48 /c-points in the irreducible Brillouin zone, 
and set Rut x ^max to 7.5, where Rut is the smallest 
muffin-tin radius and /Cmax is the plane-wave cutoff. We 
used muffin-tin radii of 2.5, 2.4 and 2.4 a.u. for Gd, Bi, 
and Pt respectively. For the non-magnetic-state calcu- 
lation valid for comparison with ARPES results at 15 
K, the seven 4/ electrons of Gd atoms were treated as 
core electrons with no net spin polarization. The atoms 



near the surface (Bi, Pt, Gd) were relaxed along the z- 
direction until the forces exerted on the atoms were less 
than 2.0 mRy/a.u.. With this optimized structure, we 
obtained self-consistency with 0.01 mRy/cell total en- 
ergy convergence. After that, we calculated the band 
structure and two dimensional Fermi surface in which we 
divided the rectangular cell connecting four ^-points by 
40 X 40, yielding 1681 /c-points. 

Even at first glance. Fig. 2 gives the impression of re- 
markable agreement between theory and experiment. All 
basic features observed by ARPES - the overall shape and 
location of the Fermi pockets [Fig. 2(a)-(b)], the binding 
energies of the bands [Fig. 2(c)-(d)] - are well reproduced 
by the calculation. The main point of this figure, how- 
ever, is the fact that band calculations show a total of 
six Fermi crossings along the f-M line segment, which 
is an even number and is not directly consistent of the 
proposed strong topological insulating phenomenoiP^I^. 
It should be noted that, in order to take into account the 
spin-orbit splitting, relativistic effects are applied to the 
calculation. Similar calculations reproduce clear topo- 
logical insulating behavior in Bi2Te3 thin films^. The 
excellent agreement shown in Fig. 2 also implies the va- 
lidity of such calculation in half-Heusler compounds. In 
fact traces for the inner two crossings is also found in 
the ARPES data, where they appear to be one single 
crossing, most likely due to finite momentum resolution 
[Leftmost part in Fig. 2(c), see also Fig. 3(d)-(h)]. 

In Fig. 3 we prove that the observed bands come from 
the sample surface. This is done by scanning the incident 
photon energy along both F-^ and f-M high symmetry 
directions. Varying the photon energy in ARPES effec- 
tively changes the momentum offset along the direction 
perpendicular to the sample surface. In our case, this 
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FIG. 4: Band structure analysis at the 
vicinity of f [red arrows in Fig. 1(e)]. 
Data is taken on LuPtBi and GdPtBi 
samples at T = 15 K. (a)-(b) Band 
dispersion maps along the f-M direc- 
tion. Green arrows point to the posi- 
tion of the inner hole band which have 
lower intensity than the two other hole 
bands, (c)-(d) Corresponding MDCs 
for panels (a) and (b). (e)-(f) Extrac- 
tion of the band position for panels (a) 
and (b). By linearly extrapolating the 
bands above the chemical potential fi 
we show an approximate band crossing 
point (Dirac point) at E ^ 0.4 eV for 
GdPtBi. 
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FIG. 5: Band structure analysis at the vicinity of M [red 
box in Fig. [ijc)]. Data is taken on LuPtBi samples, (a)-(d) 
Binding energy dependence of band structure near M. Map 
location in the surface Brillouin zone is shown in Panel (e). (f) 
Theoretical band map at the chemical potential for GdPtBi. 
(g),(h) Band maps for two perpendicular directions marked 
by red lines in (g). There are in total two Fermi crossings 
along the F-M line segment at the vicinity of M. 



direction corresponds to or the (111) direction of the 
fee Brillouin zone. Figs. 3(a)-(b) show that all observed 
bands form straight lines along the direction, a clear 
indication for the lack of kz dependence. In Fig. 3(c) we 



compare this to a calculated Fermi surface map for the 
hulk bands, along the same direction as in Fig. 3(a). The 
difference is clear: the bulk bands are dispersive along the 
V-A direction; and most of the experimentally observed 
bands are not present in the calculation. In Figs. 3(d)- 
(h) we pay special attention to the bands crossing \i near 
f by showing the band structure for four different photon 
energies. In total there are at least three Fermi contours 
surrounding F, the outer two being a lot brighter than the 
inner one (or two, see discussion for Fig. 2). As shown in 
Fig. 3(h), These three (or four) bands cross \i at exactly 
the same k positions for all photon energies. Therefore 
all of them are surface bands. The data in Fig. 3 thus 
show, unambiguously, that a metallic surface electronic 
state exists in the half-Heusler compounds. 

The exact number of Fermi crossings along the f-M 
line segment is also examined in Fig. 4. The main conclu- 
sion for Fig. 4 is that there are also three (or four) visible 
Fermi crossings at the vicinity of f between these two 
Kramers' points. We show these bands on the LuPtBi 
and GdPtBi samples. Both on the band dispersion maps 
[Figs. 4(a)-(b)] and the momentum distribution curves 
[MDCs, Figs. 4(c)- (d)] we see that there are two bright 
hole-like bands almost parallel to each other, and a much 
weaker inner band with lower Fermi velocity. This in- 
ner band is not easy to see in the band maps (nonthe- 
less indicated by green arrows), but is clearly visible in 
the MDCs by small intensity peaks tracing down from 
the one marked by a green bar [also marked by a green 
color in Figs. 4(e)-(f)]. The same band also exists in the 
V-K direction [Figs. 3(d)-(h)]. Same as the discussion 
for Figs. 2 and 3, this inner crossing is reproduced in 
the band calculation by two closely located spin-orbit- 
splitting bands that form a Kramers' pair. The brighter 
parallel bands form a second Kramers' pair of opposite 
spins. In Fig. 4(e)-(f) we show the linear extrapolation 
of the two brighter bands. In GdPtBi they are likely to 
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reduce to a Dirac point at about 0.4 eV above /i. If the 
total number of crossing is four, such a configuration wih 
give zero contribution to the total Berry phase. 

In Fig. 5 we examine the bands near the M point. The 
/c-space location of the ARPES maps [Figs. 5(a)-(d)] is 
shown in Fig. 5(e). Panels 5(g)-(h) present the band dis- 
persion maps for two cuts crossing M, whose positions 
are marked in Panel 5(f) with the band calculation result. 
Figs. 5(a)-(d) show that the M bands form a very spe- 
cial shape. At high binding energies [E ~ —0.1 eV, Fig. 
5(d)], two U-shape bands are well separated. As binding 
energy decreases these two bands merge into each other 
and hybridize to form a central elliptical contour and two 
curly-bracket-like segments. The segments near each M 
points link together, forming another large Fermi contour 
enclosing the zone center f . It is clear from Fig. 5(g)-(h) 
that there are two Fermi crossings in both the V-K and 
f-M directions. The special shape of the Fermi surface 
is formed by two bands that are likely to be members of 
another Kramers' pair. Kramers' degeneracy of spin hap- 
pens at ~ 30 meV below ji. All this features are obtained 
with our calculation for the surface states [Fig 2(b) and 
2(e)]. These two bands also give zero contribution to the 
total Berry phase. 

In summary, we performed an ARPES survey on 
the electronic structure of three half-Heusler compounds 
i?PtBi {R = Lu, Dy, Gd) which are proposed to be topo- 



logical insulators. Our result show unambiguously that 
these materials have a metallic surface state markedly 
different from the calculational result on the bulk elec- 
tronic structures. This surface state is reproduced with 
high accuracy in our band calculations. Both experiment 
and theory reveal several bands that cross the Fermi level. 
Knowledge of the exact number of these bands is possi- 
bly limited by experimental momentum resolution. No 
direct consistency with the proposed strong topological 
insulating behavior is found in the ARPES results. For 
final determination of their topological classes, both an 
APRES measurement of ultrahigh /c-resolution and a di- 
rect calculation of the first Chern number as a topological 
invariant are in need. 
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